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r~| , Recently, it has been shown that the concept of the pole mass of a 

heavy quark becomes ambiguous beyond perturbation theory, because of 
the presence of infrared renormalons. We argue that the predictions of 
heavy quark effective theory, whose construction is based on the pole mass, 
C^ ' are free of such ambiguities. In the l/mq expansion of physical quantities, 

infrared and ultraviolet renormalons compensate each other between coef- 
ficient functions and matrix elements. We trace the appearance of these 
compensations for current-induced exclusive heavy-to-heavy and heavy-to- 
light transitions, and for inclusive decays of heavy hadrons. In particular, 
we show that the structure of the heavy quark expansion is not obscured 
by renormalons, and none of the predictions of heavy quark effective theory 
are invalidated. 
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1 Introduction 

For heavy quarks, the non-relativistic bound-state picture suggests the notion of 
the pole mass rriQ"^ defined, order by order in perturbation theory, as the posi- 
tion of the singularity in the renormalized quark propagator. The pole mass is 
gauge invariant, infrared finite, and renormalization-scheme independent [IJ. In 
the context of perturbation theory, it is thus a meaningful "physical" parame- 
ter. Once non-perturbative effects are taken into account, however, this concept 
needs to be generalized, since in reality there is no pole in the quark propaga- 
tor because of confinement. Recently, it has been shown that signals for such 
non-perturbative effects can be found in the asymptotic behaviour of pertur- 
bation theory itself. The presence of infrared renormalons in the perturbative 
series that relates the pole mass to a mass defined at short distances leads to an 
unavoidable ambiguity of order Aqcd in the definition of rriQ ^ 0, ^ . The appear- 
ance of renormalons signals that perturbation theory is incomplete without the 
inclusion of non-perturbative corrections. In fact, much of the non-perturbative 
structure of a theory can be inferred from a study of the singularities of correla- 
tion functions after Borel transformation with respect to the coupling constant. 
The application of this approach to QCD was pioneered by 't Hooft [^. The 
positions of the singularities on the positive real axis signal the magnitude of 
non-perturbative corrections. In turn, the structure of non-perturbative correc- 
tions implies constraints for the structure of infrared renormalons |§]-[0- 

The existence of a "physical" definition of the mass of a heavy quark, which 
agrees with the pole mass up to terms of order Aqcd, plays a crucial role in 



the construction of the heavy quark effective theory (HQET) |p-||T^, which by 
now has become the main theoretical tool used to analyze the properties and 
decays of hadrons containing a heavy quark. In view of the intrinsic ambiguity 
in the definition of the pole mass, the question arises whether the HQET is an 
inconsistent effective theory, whose predictions are plagued by renormalon am- 
biguities. The purpose of this paper is to demonstrate that this is not the case. 
Renormalons enter HQET predictions because one tries to separate perturba- 
tive and non-perturbative (as opposed to short- and long-distance) effects into 
coefficient functions and matrix elements. Infrared renormalons appear in the 
coefficient functions since soft loop momenta give a non-negligible contribution 
to the Feynman integrals which appear in their calculation. Similarly, ultraviolet 
renormalons appear in the matrix elements because of the power divergence of 
Feynman integrals in the HQET. In this paper, we argue that in predictions for 
physical quantities such as weak decay amplitudes, renormalon ambiguities cancel 
between coefficient functions and matrix elements.[] A generic HQET prediction 

-•^For the case of the heavy quark two-point function this cancellation has been demonstrated 
in Ref. 0. 



for a physical quantity A{mQ) is of the form 

A{mQ) = Co{mQ/i2) MM + Ciimg/fi) M^ + . . . . (1) 

We will trace the cancellation of renormalons explicitly to subleading order in 
l/mq, by showing that the infrared renormalon in Co cancels against an ultra- 
violet renormalon in the matrix element Mi, so that the sum of the two terms 
on the right-hand side is unambiguous. In more complicated processes such as 
flavour-changing transitions between two heavy hadrons of different velocity, the 
way in which these cancellations take place is rather non-trivial. However, that 
they take place should not be a surprise. In fact, the appearance of renormalons 
could be avoided if in the construction of the HQET one would follow the idea 
of Wilson's operator product expansion (OPE) [l^ literally [@, |^. The OPE is 



not designed to separate perturbative and non-perturbative effects, but to dis- 
entangle the physics on different distance scales. This is not accomplished when 
one uses dimensional regularization in the calculation of the coefficient functions. 
Instead, one should introduce a hard factorization scale /i < mq by cutting out 
momenta k < ^ from the Feynman diagrams which determine the Wilson coef- 
ficients, and attribute these contributions to the matrix elements. In practice, 
this procedure is impracticable and awkward, but it would eliminate the infrared 
renormalons from the coefficient functions and the ultraviolet renormalons from 
the matrix elements. What is important is that in the HQET such a program 
could be implemented without changing the transformation property of the ef- 
fective Lagrangian under the spin-flavour symmetry [|l9l . Hence, the structure of 
the predictions obtained using the HQET remains unaffected. This implies that 
renormalons enter the usual (practical) form of the HQET in such a way that 
they do not spoil the relations imposed by heavy quark symmetry PU| and the 



equation of motion (including the vanishing of certain l/mg corrections at zero 
recoil 13), and they do not increase the number of hadronic form factors that 



appear in a given order of the l/mq expansion. 

Let us note in passing that in the lattice formulation of the HQET (or indeed 
in any regularization scheme with a dimensionful cut-off) one encounters ultra- 
violet divergences which behave as powers of the ultraviolet cutoff (i.e. inverse 
powers of the lattice spacing). These power divergences are due to the mixing of 
higher dimensional operators with lower dimensional ones. They become more 
severe as higher-order terms in the l/mq expansion are calculated. The presence 
of power divergences, and the fact that they are likely to imply the existence 
of non-perturbative effects, and hence to require non-perturbative subtractions. 



was explained in Ref . |^ . The close connection between the presence of power 
divergences and that of ultraviolet renormalons in matrix elements of higher di- 
mensional operators in the HQET was pointed out in Ref. 0, and will become 
apparent below. Techniques for the non-perturbative subtraction of the power 



divergences in lattice simulations are being developed [^]; a brief outline of the 



approach can be found in Ref. [p3 



The outline of the paper is as follows: In Sect. 0, we briefly discuss the 
appearance of renormalons in the asymptotics of perturbation theory and their 
relation to singularities in the Borel transform of correlation functions. To obtain 
a renormalon calculus which is convenient for explicit calculations, we follow 
Refs. [@, 0] and consider QCD in the limit of a large number of light quark 
flavours. We then recall some of the reasoning behind the usual construction of 
the HQET and show in which way it is affected by infrared renormalons in the 
pole mass of the heavy quark. In Sects. |^ and H, we study the cancellation of 
renormalons in exclusive heavy-to-heavy and heavy-to-light decay processes. In 
the first type of decays, the symmetries of the effective theory imply a set of 
non-trivial consistency conditions, which relate the infrared renormalons in the 
coefficient functions of bilinear heavy quark currents to the infrared renormalon in 
the pole mass. We derive the exact form of these relations, which are independent 
of any unknown hadronic matrix element. We then check them to order l/Nj. We 
also show with an explicit calculation that a sum rule recently derived by Shifman 



et al. ||2J] , which has been used to put a bound on the hadronic form factor that 
enters the extraction of \Vcb\ from semileptonic decays, cannot be correct, due 
to a mismatch of infrared and ultraviolet renormalons. In Sect. ^ we show that 
renormalon contributions cancel in inclusive, current-induced decays of hadrons 
containing a heavy quark. This proves a conjecture of Bigi et al. |@, although 
we do not agree on the details of the cancellation. In Sect. ^, we summarize our 
results and give some conclusions. 

2 Renormalons and the Construction of the HQET 

Given a perturbative series for some quantity F{as) in terms of the coupling 
constant as{^) renormalized at some scale fi, 

F{as) = T.Fn[^as{f.)j , (2) 

where /3o = H — ^Nj is the first coefficient of the /3-function, we define the Borel 
transform F{u) of F{as) by 

CO -| 

Fiu)=Fo5iu) + J2-F^+,u''. (3) 

If the series is Borel summable, the function F{as) can be reconstructed from its 
Borel transform using the integral relation 

oo 

F(„,)=/d,.exp(-^)FM. (4) 



However, if the coefficients Fn in (Q) develop a factorial divergence for large n, the 
Borel transform F{u) can have singularities on the integration contour, and the 
naive Borel summation fails. In such a case, the result of the integration depends 
on a regularization (or resummation) prescription, and F{as) is not uniquely 
defined in terms of F{u). 

In QCD, one source of divergence in the expansion coefficients of a pertur- 
bative series is related to higher-order diagrams in which a virtual gluon line 
with momentum k is dressed by a number of fermion, gluon and ghost loops. ^ 
Effectively, this introduces the running coupling constant gs{k) at the vertices. 
Since the coupling constant increases for low momenta because of asymptotic 
freedom, the insertion of additional bubbles drives the gluon line to increasingly 
softer momentum, i.e. the infrared region in Feynman integrals becomes more 
important. When the running coupling constant is expressed in terms of a fixed 
coupling constant renormalized at some large scale /i, using 

l-|a.(,)l4 - ^^^ *^ 

the appearance of powers of large logarithms leads to a factorial divergence in the 
expansion coefficients Fn in (^. Associated with this are renormalon singularities 
in the Borel transform F{u). 

In our case, the renormalon singularities will occur as single poles on the real 
axis in the Borel plane. Poles on the positive real axis, which arise from the 
low-momentum region of Feynman diagrams, are called infrared renormalons.0 
Let us denote the positions of these poles by Ui and their residues by Tj, so that 

^M=E^^ + --- (6) 

where the ellipses represent terms that are regular for m > 0. For the calculation 
of the inverse Borel transform from (§) , we may write the pole denominators in 
terms of a principle value and a (5-function contribution: 

1 1 

^ P + Vi Ku - Ui) . (7) 

U — Ui U — Ui 

Here, rji is a complex number which depends on the regularization prescription. 
For instance, one may choose one of the following regularizations (with 5 -^ +0) 



u — Ui + k6 
'u — UiY + 6'^ 



Vi 



^We hasten to add that in a non-abehan theory bubble summation is not a gauge-invariant 
procedure. This is one of the reasons why we will have to use a large- A^/ expansion to obtain 
a consistent renormalon calculus, see below. 

^Similarly, poles on the negative real axis arise from the high- momentum region and are 
called ultraviolet renormalons. 



1 

— ^ -^ Vi = ±«7r . (8) 

U — Ui^ to 

One may also choose the principal value prescription itself, in which case rji = 0. 
We write the regularized form of the Borel transform as 

F{u) = F.egiu) +J2'niri5iu-Ui) , (9) 

i 

where by definition Fj.eg{u) contains the pole terms regularized with a principle 
value prescription, and all ambiguity resulting from the freedom to use a different 
prescription resides in the ^-function contributions. The inverse Borel transfor- 
mation then leads to 

F{as) = Frcg(as) + ^Viri exp ( - -^ ) ~ Freg(a,) +^Vir^( -^^ 

(10) 
where -Freg(as) is the inverse Borel transform of Fj-eg^u). In the last step, we have 
used the one-loop expression 

"■<"' = /?„ ln(,VA'5co) '''' 

for the running coupling constant .0 These definitions make explicit the fact that 
terms which depend on the regularization prescription are exponentially small in 
the coupling constant, i.e. they have the form of power corrections. The leading 
asymptotic behaviour is determined by the nearest infrared renormalon pole at 
u = ui. We define the renormalon ambiguity AF as the coefficient of rji. 



AF = n exp - ^ ~ n ^^ . (12) 

It is a measure of the intrinsic ambiguity in the quantity F arising from the 
necessity to regularize the divergent behaviour of perturbation theory in large 
orders. It is the purpose of this paper to trace how these leading (in powers of 
Aqcd//^) ambiguities cancel in HQET predictions for physical quantities. 

Although the appearance of renormalons can hardly be doubted on physical 
grounds, a rigorous proof of their existence does not exist even in field theories 
that are much simpler than QCD. For this reason, various forms of large- A^ expan- 
sions have become the state-of-the-art approach to study renormalon singulari- 
ties. In QCD, one uses 1/A^/ as an expansion parameter, where Nf is the number 
of light quark fiavours. In the large- iV/ limit, the insertions of fermion loops in a 
gluon propagator are the only higher-order contributions that have to be retained 
in the perturbative expansion, since they involve powers of Nf a^ = 0{Nj). Un- 
fortunately, QCD in the large- iVj limit is not an asymptotically free theory; the 

''Note that the last relations in (O) and (ll^) become exact in the \arge-Nf limit. 



first coefficient of tlie /3-function becomes negative for Nf > 33/2. However, it 
is believed that although the 1/A^/ expansion is not adequate to describe the 
dynamics of QCD, it can still be used to locate the position of the renormalon 
poles in the Borel plane. In other words, the hope is that tracing the fermionic 
contribution to the /9-function one gets the remaining contributions for free, and 
that using the correct value of /3o in (§) gives the right result. Although there 
exists no proof of this assertion, we will accept it as a working hypothesis. 

In the large- A^j limit, the summation of bubbles can be performed directly on 
the gluon propagator. In Landau gauge, and after renormalization of the fermion 
loops, the Borel transform of the resummed propagator takes the form 0, |^ 

where /i is the renormalization scale, and C is a scheme- dependent constant. 
In the MS scheme, C = —5/3. Consider now an arbitrary correlation function 
without external gluons. To order l/Nf, all its dependence on the coupling 
constant a^ comes from diagrams containing one resummed gluon propagator. 
The Borel transform of such diagrams is simply obtained by using the Borel 
transformed propagator (|T3|) instead of the usual propagator. 

Following Beneke and Braun , let us then consider the structure of infrared 
renormalons in the pole mass and on-shell wave-function renormalization of a 
heavy quark. In terms of the self-energy S(^, one has 



pole 



^Q =^Q + S( 



-1 d^lD 



1 ^Q -^ 



(14) 



where mq is the bare mass. In general, these are complicated implicit equations. 
However, since the self-energy is of order 1/iVj with respect to the bare mass, 
one can replace m^ ^ by rriQ on the right-hand side, thereby neglecting terms of 
order 1/A^/- We work in Landau gauge and use dimensional regularization. By 
evaluating the diagram depicted in Fig. |1], we obtain for the Borel transform of 
the self-energy the relations 



E(V, u) L = ^ [d-l) e-^" [A.r'^l^ (^) 



d-A-2u 



x(,-2-2. /'^-it,"'^"'-f-^"' .0(iVr). 



9S(j/,«) 



d^ 






T{d-l-u) 
S(j/,n), +0{N7'), (15) 



=mQ 



where Cp = 4/3, and d denotes the number of space-time dimensions. For generic 
u, one can evaluate these expressions for d = 4. The positions of renormalons are 



determined by the F-functions in the numerator. There are infrared renormalons 
at positive half-integer values of u, as well as ultraviolet renormalons at negative 
integer values of u. For u = 0, the self-energy and its derivative are ultraviolet 
divergent in d = 4 dimensions. One can subtract the ultraviolet divergence by 
subtracting the pole at m = after setting d = 4 0. This determines the 
renormalized Borel transform up to a scheme- dependent function R{u), which is 
entire in the Borel plane if a renormalization scheme with analytic counterterms 
(such as MS) is employed. The result is 



— pole / \ 



Ulr 



Sin) + ^ 



+ 0{Nf 






2u 



U 



T{u)T{l-2u) 



F(3 



u 



3 

u 



- + Rm{u) 



(16) 



and 



JR, 



u 



Siu) + ^ 






2m 



U 



nu)ni-2u) 3 

F(3-n) +^ + ^^(^) 



+ 0{Ny^) 



(17) 



where rriQ is the renormalized mass. The first expression has been derived in 
Ref. 10]. Note that to order 1/Nf the choice of tuq in the parentheses on the 
right-hand side is arbitrary. The functions Rmiu) and Rziu) depend on the 
renormalization scheme specified by the superscript "R" . In the MS scheme, one 
has Rmiu) = —5/2 + 0{u) and Rziu) = 11/2 -|- 0{u). For the discussion of 
renormalon singularities these functions are irrelevant. The asymptotic behaviour 
of the perturbative expansions for Tng'^ and Zq is determined by the nearest 
infrared renormalon pole, which is located at m = 1/2. According to (|T2[), it 
leads to intrinsic ambiguities given by 



pole 



AZ 



_2Cf 
3Cf 



e-^'/^ Aqcd + 0{Ny 



Ac 



Po mq 



f) 



Note that the product e~'"/^AQCD is scheme-independent. 

After this lengthy introduction into the problem, let us now turn to the con- 
struction of the HQET |^-[|ll|]. A heavy quark interacting with light degrees of 
freedom inside a hadron is almost on-shell. It is then natural to split its mo- 
mentum into a "large" and a "small" piece according to pq = mq v + k, where 
V is the velocity of the hadron, and rriq is some choice of the heavy quark mass 
discussed in detail below. For the moment let us just require that the components 
of the residual momentum k are much smaller than mq. One then proceeds by 
introducing a velocity-dependent heavy quark field Kix), which is related to the 



original field Q{x) by 



(1 +^) 
hy{x) = exp{imQ v ■ x) — - — Q{x) . (19) 



The effective Lagrangian for hy reads [T^, |Tl|, [15 



-^cff = hy {iv ■ D — 5m) hy + . . . , (20) 

where 6m is the residual mass term for the heavy quark in the effective theory. 
It appears since there is a freedom in the choice of the expansion parameter 
m,Q in (|r^) . One can show that in physical matrix elements only the combination 
{m,Q+6m,) appears, i.e. different choices oimg are compensated by different values 
of 6m, ||15|. The ellipses in (pOf) represent terms that contain additional powers of 



iD^/m,Q or 6m,/ mq. If one arranges things in such a way that the components of 
k and 6m are of order Aqcd and independent of m,Q, this construction provides 
a systematic expansion in powers of Aqcd /m,Q. Moreover, the leading terms in 
the effective Lagrangian (^) are then invariant under a spin-flavour symmetry 
group. To this end, the heavy quark mass m,Q used in the field redefinition 
( P^D must be a "physical" mass such as the pole mass, the mass of the lightest 
hadron that contains the heavy quark, or any other definition that differs from 
the pole mass by an amount of order Aqcd- The residual mass term is given 
by 6m, = itiq"^ — m,Q; i.e. if one chooses the pole mass to construct the HQET, 
the residual mass vanishes, and to any finite order in perturbation theory the 
effective heavy quark propagator has a pole at k = 0. However, from our previous 
considerations we know that there is an intrinsic ambiguity of order Aqcd in the 
definition of the pole mass, once non-perturbative effects are taken into account. 
Hence, if one wants to write down the Lagrangian of the HQET without specifying 
a particular Borel summation prescription, one can do this for the price of an 
ambiguous residual mass term 0]. To be specific, let us construct the HQET using 
the heavy quark mass defined with a principle value prescription to regularize the 
poles in the Borel plane [cf. (^^]. It then follows that 

6m = vi AmP°^^ = -r^i ^ e'^/^ Aqcd + 0{Nj') . (21) 

Po 

The ambiguity associated with the definition of the pole mass shows up in the 
form of an ambiguous parameter in the effective Lagrangian (pOl) . At first sight 
this may seem a problem: How can one derive unambiguous predictions from 
a Lagrangian that contains an ambiguous parameter? The answer is that the 
effective theory has to be matched onto the full theory at some large momentum 
scale. In this process there appear coefficient functions multiplying the operators 
of the HQET. The ambiguous residual mass term is required to cancel ambigu- 
ities in these coefficient functions. The important point to note is that 6mL is 
independent of ulq and thus does not break the flavour symmetry of the effective 

8 



Lagrangian. The way in which the residual mass enters the l/mg expansion has 



been investigated in Ref. |T5 



Let us come back, at this point, to the original formulation of Wilson's OPE 
1^ , in which renormalons never appear. Introducing a hard factorization scale 
yU in the construction of the HQET would yield a residual mass term of the form 
6171 ~ fias{fi). Likewise, hadronic matrix elements in the effective theory as well 
as the Wilson coefficient functions would have a power-like dependence on /i, in 
such a way that the factorization scale disappears from the final predictions for 
physical quantities. This is the content of the renormalization-group equation. 
In this formulation, the parameters of the theory are not plagued by ambiguities, 
but they depend on the arbitrary parameter /i. Moreover, the precise form of 
this dependence (for instance, the coefficient of the fiasin) term in 6m) depends 
on how exactly the hard cutoff is implemented in Feynman diagrams. Hence, 
there is a similar arbitrariness in the definition of these parameters as in the case 
of the practical form of the OPE, which contains renormalons. Finally, we note 
that chosing /i = rriQ as a factorization scale (as it was proposed in Refs. 0, p3| ) 
breaks the fiavour symmetry of the effective Lagrangian (through the residual 
mass term) and is thus not a viable choice in processes that involve more than 
one heavy quark fiavour. 

3 Heavy-to-Heavy Transition Matrix Elements 

In this section we investigate how renormalons appear in the hadronic matrix ele- 
ments that describe current-induced transitions between two hadrons containing 
heavy quarks with masses mi and m,2- These matrix elements are of the form 
{H2{v2)\ Q2^ Qi \Hi{vi)), where vi and V2 denote the velocities of the hadrons. 
The quantum numbers of the light degrees of freedom are assumed to be the 
same in the initial and final state, but are otherwise arbitrary. We will restrict 
ourselves to the cases of vector or axial vector currents (F = 7^^ or 7'^75). This 
covers semileptonic weak decays such as B ^ D^*^^ and A^ ^ Aciu. 

In the HQET, the currents which mediate these transitions obey an expansion 
in a series of local operators multiplied by coefficient functions. These functions 
depend upon the heavy quark masses, the renormalization scale, and the hadron 
velocity product w = Vi ■ V2- A particular property of heavy-to-heavy transi- 
tions is that the coefficients of the operators of dimension four are all related 



to the coefficients of the dimension-three operators |26|. The reason for this is 



an invariance of the effective theory under reparametrization of the heavy quark 



momentum |T6|. As an example, we give the exact form of the expansion of the 



vector current to order l/mq 

Q2 i^Qi -^ CX \ Kj'^K, + -^ K,^^iTpi K, - -^ K, ip2 YK, 

I zmi 2m,2 



ow yrrii 1712 J 

+ C2 \ hv2 Vi K-, + 7^ hv2 Vi %P\ K^ - - — h^^ iTp2 Vi K^ H h^^ iV^ h^^ 

y znii 2m2 nil 

dCY i 1 J „. ^, It.^u,! 

+ -7: — i — h^2 v'{iV2-Vihy^ /i^2 ivi-V2v'[K^ ) 

ow yrrii 1712 J 



1 . „ .^ , 1 . .t^ „ . 1 



+ C3 I K^ v^ K, + - — K^ v!^ iTpi K^ - - — /i^2 iTp2 f 2 K^ V ^^2 h. 



/mi lm2 m2 



VI 



H — 

dw 



I — K^ v!^iv2-'Dihy^ /i^2 ivi-V2v!^hy^ } 

[mi m2 J 



+ 0(—„—,,^—), (22) 

where CY = CY{mi/fi,m2/fi,w). A similar expansion with coefficients Cf, and 
with 75 inserted after whatever object carries the Lorentz index /i, holds for the 
axial vector current. The symbols T>i represent combinations of a gauge-covariant 
derivative and the residual mass term. They are defined as [^ 



iV^ = iD^ -57nv^, iV^ = iD^" + 67nv^. (23) 

We will show below that the coefficient functions Q ' of the dimension-three 
operators contain infrared renormalons at m = 1/2, corresponding to power be- 
haviour of order l/rrii or l/m2- In order for the physical heavy-to-heavy transi- 
tion amplitudes to be unambiguous, we have to require that these renormalons be 
compensated by ultraviolet renormalons in HQET matrix elements of dimension- 
four operators. This requirement is analogous to the renormalization-group equa- 
tion in Wilson's OPE. The complete set of dimension-four operators consists of 
the local current operators in ( p2l) as well as operators containing the time-ordered 
product of a dimension-three operator with a I/tuq insertion from the effective 
Lagrangian |]14|, which at this order is given by |]8|, [TT| 



jCe{i = Kiv-Vhy + - hy{iVfhy + Cmi,gi'fnQ/n)—^—hya^uG^''K + ... . (24) 

However, a cancellation of renormalon ambiguities can only occur between terms 
that have the structure of matrix elements of local dimension-three operators. In 
other words, only the matrix elements of dimension-four operators that can mix 
with lower dimensional operators can contain ultraviolet renormalons. 

In heavy-to-heavy transitions, the ultraviolet renormalons in the matrix ele- 
ments of the local dimension-four operators can be related to the infrared renor- 
malon in the pole mass. Using the equation of motion of the HQET, iviT>i hy^ = 0, 
as well as an integration by parts, one can show that 0, ^ 

{H2{V2)\ K.TiV'iK, \Hi{vi)) = ^^ {H2{v2)\ K, r (w< - t;°) K, \Hi{vi)) + . . . , 

10 



-{H,{v2)\ K,t^rK^\H,{v,)) = J^{H2iv,)\ K,{wv^ - o r/i^ji/iK)) + . . . 



(25) 



where F may be an arbitrary Dirac matrix. The eUipses represent terms that 
vanish upon contraction with Via or V2a- These terms cannot be written in the 
form of matrix elements of local operators. Hence, as explained above, they must 
be free of renormalons. As an example, consider the case of the ground-state 
pseudoscalar and vector mesons. There, the matrix elements of local dimension- 
three operators can be parametrized in terms of the Isgur- Wise function pOf : 

{M2{V2)\KTK, \M,{v,)) = -aw,fi)Tr{M2{v2)TMi{vi)}, (26) 

where A4(v) are the tensor wave functions defined in Ref. [O. The matrix 



elements of local dimension-four operators can be written as 0, ^ 

{M2{v2)\K,rtV'^K,\M^{v,)) = ^aw,f^)TT{M2{v2){wv'^-v^)TMl{v^)} 

w + 1 '- J 

+ e3(^, /i) Trj (7" - ^£±^) 7W2(.;2) r A<i(t;i)} • 

(27) 

Note that the Feynman rules of the HQET imply that there cannot appear Dirac 
matrices next to F under the trace with the meson wave functions. Obviously, the 
structure of the trace associated with the function ^3(1/7, /x) is different from the 
structure of the trace in (|26|) . It follows that C.si'w, fi) does not contain an ultravi- 
olet renormalon at m = 1/2. Let us now come back to the terms shown explicitly 
in (|25|). They have the structure of matrix elements of the local dimension-three 
operators. For instance, in the case of the second operator on the right-hand side 
in (^2|) one has F = 7'^7a, and between the heavy quark spinors one can replace 
{wvi — V2) 7'^ 7a by {w + l)'~f^ — 2^2. The parameter 

A = rriHi — rrii — 5m = rriH. — m\° ^ ; i = 1,2 (28) 

denotes the asymptotic value of the difference between the hadron and heavy 
quark pole masses, which is flavour-independent. Note that this parameter is 
independent of the choice of the expansion parameter rriQ used in the construction 
of the HQET [|l^. However, because of its dependence on the pole mass it does 



contain an ultraviolet renormalon [0. Using (plf), we find that the corresponding 
ambiguity in A is given by 

AA = -AmP°'^ = -^ e-^/2 Aqcd + 0{Ny') . (29) 

Next consider the matrix elements of the operators containing the time- 
ordered product of two local operators. An insertion of the kinetic operator 
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(l/2mQ) hy {iT>Yhv into a matrix element of a local dimension-three operator 
does not affect the transformation properties under the Lorentz group and heavy 
quark spin symmetry. The effect of such an insertion is simply a multiplicative 
renormalization of the original matrix element. We define a function K{w, fi) byQ 

= K{w,fi) {H2iv2)\ /i,,(0)r/i,,(0) \H,{v,)) . (30) 

Clearly, these time-ordered products can mix with the dimension-three operators 
under renormalization. This is obvious if a dimensionful regulator is employed. 
But even in dimensional regularization, it can be seen from the renormalization- 
group equation for the function K{w, /x), which contains an inhomogeneous term 
proportional to A P7| : 

/i -^ K{w, /i) = -2A («; - 1) ^ 7hh(«^) . (31) 

d/i aw 

Here 7hh(w^) denotes the velocity-dependent anomalous dimension of the currents 
in the effective theory |T^: 



Thhlw^J 



wr{w) — 1 + 0(a^ 



vr 
1 



sJ 1 



r{w) = ^ In {w + Vw2 - 1) . (32) 

V w^ — 1 

From (0), it follows that the function K{w, fi) contains an ultraviolet renormalon 
at M = 1/2. Let us denote its renormalon ambiguity by AK{w). Vector current 
conservation implies that K{w,fi) must vanish at zero recoil {w = 1) |]l4l, and 
this requires that 

AK{1)=0. (33) 

Finally, we note that insertions of the chromo-magnetic operator change the 
transformation properties of matrix elements in such a way that there is no 



mixing with matrix elements of lower dimensional operators [|I4|. Hence, the 
HQET functions that parameterize these matrix elements are free of ultraviolet 
renormalons. This is again special to the case of heavy-to-heavy transitions, 
where the spin symmetry applies to both the initial and final state. 

We can now equate the infrared and ultraviolet renormalon ambiguities to 
derive a set of conditions that have to be fulfilled in order to obtain unambiguous 
predictions for the physical heavy-to-heavy transition form factors. Separating 



^For meson decays, the form factor K(w, /i) is usually written as K(w, /i) 
2xi{w,^)/^{w,fi) iQ, where f(w,//) is the Isgur-Wise function pC| ]. 
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the terms associated with different Lorentz structures, we obtain from (^2j), (p5D 
191) and (|3^ 



1 \iTii 1122/ lw + 1 OW ) 2 Krrii m2, 

AC, = ^^^ f^ + ^V^^^C, + 2(u; - 1) ^1 + ^^^ ^^^C 

2 Vmi m2/lw + l aw J mi w + 1 



?7l2 W + 1 2 Vmi 7712/ 

2 Vmi m2/lw + l ow ) m2 w + 1 

T^^4.(C.±C. + C3)-^(± + ±)c,. (34) 

mi w + 1 2 Vmi m2/ 

We use a short-hand notation where we omit the superscript V or A on the 
coefficient functions, and where upper (lower) signs refers to the coefficients of 
the vector (axial vector) current. In total, there are thus six relations. We find 
it useful to solve for AK using the relation for AC^, and to eliminate AK from 
the remaining relations. This leads to 

Air = Ampoie(^-l)^Y + 2^1nCi^ -2 — + — — ^ , (35) 

lw + 1 ow J Vmi m2/ Cf 



as well as 

1 (ACX 
Ampoie V CY 

' (AC- 

Ampoie V 

/ (AC,- 

Ampoie V 


AC^ 
Cf J 


-rS M^ ^ 1 (w I) ^ \n^^^ 


Ui ' m2JU + l ' ^"^ ^hw^'^Cfi' 

- ' i rci±C2 + C3)+ ^ "^1^2 

m2 w + 1 mi U7 + 1 
Vmi m2/ ow Ci 

mi w + I in2 w + L 

+ (' + M(- 1)^3^ ln^^ (36) 
Vmi m2/ ow Ci 



Equation ( pSf ) determines the structure of the ultraviolet renormalon in the 
hadronic form factor K{w, fi) in terms of the infrared renormalon in the pole 
mass and in the coefficient function Cj^. Since we are not able to calculate the 
hadronic form factor K{w, n) from first principles (not even using a 1/A^/ expan- 
sion), we cannot check this relation, but we can use it to compute AK given a 
calculation of Ampoie and ACj^. The remaining five relations in (^) form a set of 
consistency conditions involving only the infrared renormalons in the pole mass 
and in the coefficient functions. Unless these conditions are satisfied, a compen- 
sation of infrared and ultraviolet renormalons is not possible. The existence of 
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such relations is non-trivial and is a consequence of the strong constraints im- 
posed by heavy quark symmetry on the structure of the weak decay form factors 
in heavy-to-heavy transitions. 

The above results are exact to all orders in 1/Nf. In the large- A^/ limit, they 
simplify since the renormalon ambiguities are of order l/Nj, and we can use the 
fact that Ci = ^n + 0{1/Nf). This leads to 



AK 



w — 1 



Am 



w 



pole 



— + — ) AC^ + OiNf 
nil ^12/ 



(37) 



and 



ACX - AC^ 



Am 



pole 



W 



V.A 



AC. 



ACV 



1 



- + -) 

mil fn2/ 
Ampoi, 



0{N~r 



T 



w + 1 m,2 

1 Ampole 



w 



1 mi 



+ 0{NJ') 
+ 0{Nf) 



(38) 



Let us now check these relations with an explicit calculation of the asymptotic 
behaviour of the coefficient functions to order l/Nj. The coefficients are obtained 
by comparing matrix elements in the HQET with matrix elements in the full 
theory at some reference scale fi. This matching procedure is independent of the 
external states, and it is most economic to evaluate the matrix elements with 
on-shell quark states. If one uses dimensional regularization, all loop diagrams in 
the HQET vanish, i.e. the coefficient functions are simply given by the on-shell 
vertex functions of the full theory 
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Hence, the Borel-transformed coefficient 
functions are obtained by evaluating the diagram shown in Fig. || supplemented 
by wave- function renormalization. Setting d = 4 in the final result, we obtain 



W,A 



C^'^'iu) = 5iu) 



3C 






F^-Cu 



2u / II \ 2«1 



+ 



/i- 



— 

m^iJ \mi2J 

\"r(M)r(i-2M) 



(1 



u 



mim,2 



J r(2 



u 



r(3 - u) 

;1 + m)w;±m1Fi\+" 



Co (u) 



c^ 

_2Cf 



+ u 
2 r 



^1 Fl+« 
m2 



^2 pi+„ 
mi 



+ 2 



{l-u){l-2u) 
(2-n) 



Fn 



u 



w r{w) - 1 + RvA^) \ + 0{Nf) 



-Cu 



/ /i^ \"r(i + n)r(i-2M) 



Fi+« 



\mt1m2J 
{l-2u 



3(2 



u 



22 



r(2-'u) 

2mLi 



± 



C3' (m) 



T^e-^" 



/5o 



/i^ 



m,im2/ 



m2 

yr(i + M)r(i 



pl + M 



+ 0{NJ 



2u) 



r(2 



u 
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Fr - 1^ (i^r ± ^ ^-" ) K ^^^") ' ^39) 




where 



1 



F„^ -Ii^±^ /dxx^-Hl-a;)''- 



afe 







a;2!!^ + (l_a;)2!!!a + 2x(l-x)w 



(40) 



r(a)r(6) 

ms in the 
olet divergences for -u = 0. The coefficient of the 1/u pole is proportional to the 



The terms in the last row in C^ ' {u) come from a renormalization of the ultravi 



one-loop coefficient of the velocity-dependent anomalous dimension in (32). The 
detailed form of the entire function i?y_^(-u) is irrelevant for our discussion. We 
note that a check of the complicated expressions (|39|) is provided by an expan- 
sion around m = 0, from which we recover the one-loop results for the coefficient 
functions derived in Ref. p8| . 

It is a simple exercise to extract the residues of the renormalon poles at 
u = 1/2 from the above expressions. The relevant parameter integrals are given 
by 

FT = ^ ./^, FT = 4t ^/^' ^-' = \ (^- VF,f ) . (41) 
w + 1 \} 1712 w + 1 y mi 2 




Using (IT4) and (|T8D, we then compute the corresponding renormalon ambiguities 
in the coefficient functions. We find that the consistency conditions (|3^) are 
indeed satisfied. In particular, we note that 

ACr = ^^^ ^^^^ (- + -)+ 0{Nf) , 
4 w + 1 Vmi ma/ ^ / ^' 



^^ _ Ampole 3{W-1) / I 1 \ ^,,.-2x 



^ 4 w + 1 Vmi 7712 / ^ ^ ^' ^ ^ 



so that the difference obeys the ffist relation given in (^Sj). Moreover, for the 
ultraviolet renormalon ambiguity in the function K[w, n) we obtain from 

w -1 Amp ole , r^rAT-2\ 
'w + 1 2 



AK{w) = -^^ =3^ + OiNy') . (43) 



Note that AK{w) vanishes for w = 1, as required by vector current conservation 
[see (ID]. 

We have emphasized above that the appearance of renormalons does not ob- 
scure the structure of the heavy quark expansion. In particular, Luke's theorem 



IJ] , which concerns the vanishing of ffist-order power corrections at zero recoil, 
remains unaffected. Let us illustrate this important fact with the example of the 
meson decay form factor hAi{w) defined by |jlj 



{D*{v2, e)| C7^75& \B{vi)) = ^mBmo* [w + 1) hA^w) e*^ + . . . , (44) 
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where w = vi ■ V2-, and e denotes the polarization vector of the Z)*-nieson. This 
form factor plays a crucial role in the extraction of | Vcb\ from the extrapolation 



of the B -^ D*£h' decay rate to zero recoil. One obtains [22 



1 dT(B^D*iu) Gl , ,, o , ,o,, , ,,o 

hm -^== -^—. ^ = -^ ms - mo^fml, V,,' /^^, 1 ^ . 45 

"'^i y/w^ — 1 dt(7 47r'^ 

The important point is that /iAi(l) is protected by Luke's theorem against first- 
order power corrections [Q. It follows that 



hA,{\) = r/^ + 0(l/my ; r/^ = C^{w = 1) , (46) 

where we use rriQ as a generic notation for rric or rrih. The presence of an in- 
frared renormalon at m = 1/2 in the short- distance coefficient rjA would spoil 
this non-renormalization theorem. However, from ( ^7|) and (33) it follows that 
the infrared renormalon at m = 1/2 in C^ vanishes at zero recoil. Our explicit 
result (1^) confirms this to order 1/Nf. Thus, the theoretical uncertainty in the 
determination of \Vcb\ is of order l/wig; it is not affected by a renormalon am- 
biguity of order l/niq. Note, however, that the expression for Ci{u) in (^) 
contains a renormalon pole at m = 1, which does not vanish at zero recoil. The 
corresponding ambiguity in rj^ is given by 

Ar/^ = ^e-^A^CDf— + -V + 0(iV7^). (47) 



r,n - --QCD ^ -T^^l^f 

2po ^ \mi 1712/ 

This infrared renormalon must be compensated by an ultraviolet renormalon in 
the terms of order l/rriQ in (^Bf). For completeness, let us also study the renor- 
malization of the vector current at zero recoil. There, the relevant combination 
of coefficient functions is 

r]v = CYiw = 1) + C^iw = 1) + Cliw = 1) . (48) 

We find that the leading renormalon pole in the Borel transform of riv is located 
at M = 1. From its residue, we obtain 



^- -CM^fJ__l_\ 



2 



Note that Ar]v vanishes in the limit rrii = m2, in which the vector current is 
conserved and not renormalized, and hence r]y = 1. 

As a second example, we demonstrate the cancellation of renormalon ambi- 
guities in the ratio of the vector form factor hyiw) defined by 

{D*{v2,e)\-c-i^h\B{v{)) = ^mBmn*hviw)e^'"^^e:v2aVip (50) 
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and the axial form factor /lyi^(w). Including power corrections of order l/nic (and 
neglecting those of order l/r/i;, and higher), one finds that [^ 



hv{w) _CY{w) f^^^ 



+ («^ - 1) TTI In 



d ^cr{w 



(51) 



hA^iw) Ciiw) I rricVw + l dw Ciiw) 

Using the first relation in (0), we see that the infrared renormalon of order 1/mc 
in the ratio CY /C^ is precisely compensated by the ultraviolet renormalon of the 
term proportional to A. 

We conclude this section by pointing out an important implication of our 
result (^71). Recently, it has been claimed that one can derive a sum rule for 
the form factor h^ii^)-, from which it is possible to obtain a bound for the non- 
perturbative corrections of order I/toq in (^61). The sum rule reads p4 



2 ^2 Ai + 3A2 / 1 1 ^ K 3 



fcl(l) + ...=-<i-^ + ^l-^t,;^ + ^ + 5;^j+0(l/™J), (52) 

where the ellipses represents positive contributions from transitions into excited 
states. The parameters Ai and A2 are defined in terms of the 5-meson matrix 
elements of the kinetic and the chromo- magnetic operator in the effective La- 



grangian (|2J). In principle, these HQET parameters could contain ultraviolet 
renormalons. However, since A2 is proportional to the mass splitting between B 
and B* mesons, it is protected from renormalons. Moreover, from the expansion 
of the meson mass itib in powers of l/mi, (this extends ( pSj) to order l/nii, pO[] ) 

ms = m^ + A - h±3h + 0{l/ml) , (53) 

and from the fact that to order 1/A^/ the Borel transform of the pole mass given 
in (|T^) does not contain an infrared renormalon pole at m = 1, it follows that 
Ai does not contain an ultraviolet renormalonQ (at least) to order 1/Nf 0. We 
conclude that the non-perturbative corrections in (^) do not contain the ultravi- 
olet renormalons required to cancel the infrared renormalon in the perturbative 
coefficient r]\. Hence, there must be something wrong with the sum rule. Ei- 
ther the short- distance correction on the right-hand side is not given by r/^, or 
there must be additional terms of order l/wig to compensate the renormalon in 
ri\. Therefore, the numerical implications derived from this sum rule in Ref. pi 
should be taken with caution. 



4 Heavy-to-Light Transition Matrix Elements 

We have seen in the previous section that the appearance of renormalons in 
heavy-to-heavy transition matrix elements is to a large extent constrained by the 

^Note that even if Ai and A2 would contain ultraviolet renormalons, the mass dependence 
of the power corrections in (p2) would not match with the mass dependence of the infrared 



renormalon pole in tja as given in (47) 
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symmetries and equation of motion of the HQET, which apply to both the initial 
and final hadron states. As a consequence, ultraviolet renormalons enter the 
HQET matrix elements of dimension-four operators only through the parameter 
A and a single function K{w,fi). Since the ultraviolet renormalon in A is related 
to the infrared renormalon in the pole mass, it is possible to derive the consistency 
relations (|36|), which determine the infrared renormalon poles in the coefficient 
functions independently of any unknown hadronic matrix element. 

It is well-known that in heavy-to-light transitions there are fewer constraints 
imposed by heavy quark symmetry. In particular, most (if not all) form factors 
appearing at order l/mq mix with lower dimensional operators under renor- 
malization. Examples are provided by the l/mq expansions for meson decay 
constants [^ and the semileptonic B ^ tt ii? decay form factors |32|. Therefore, 



it is not possible to derive consistency relations analogous to (^) in this case. 
The best one can achieve is to deduce the structure of ultraviolet renormalons in 
the hadronic form factors of the HQET from a calculation of the infrared renor- 
malons in the coefficient functions and the pole mass. We shall discuss this for 
the simplest case of meson decay constants. 

Consider heavy-to-light transition matrix elements of the form {X\ qT Q\H[v)), 
where F = 7^^ or 7^75, H{v) is a heavy hadron with velocity v, and X is some 
light final state. For simplicity, we set the mass of the light quark to zero and use 
a regularization scheme with anticommuting 75. This leads to a simple relation 
between the coefficient functions appearing in the expansion of the vector and 
axial vector currents |T^: 

qYQ -^ Ciimg/fx) qYK + C2{mQ/fi) qv^K + 0{l/mQ) , 

qi^l5 Q -^ Ciimg/fi) qYl5 K - C2(mQ//i) qv^'-f^ K + (^(l/mg) . (54) 

These coefficients can be calculated in analogy to the previous section. For their 
Borel transforms, we obtain 



CAu) 




F(m)F(1-2m) 
F(3-m) 



r(\ 4^^ -cnf /^ V^ T{l + u)T{l-2u) , ^,^^_,, 
Po \mQj F(3-m) ^ 

where R{u) = 5/4 + 0{u) in the MS scheme. We have checked that from an 
expansion around u = one recovers the known one-loop expressions for the 



coefficient functions given in Ref. |^. From the residues of the poles at u = 1/2, 



it is straightforward to compute the renormalon ambiguities 

12 mq ^ 



2 Ampoie 

3 niQ 



AC2 = ^ =^^^^ + 0{Nf) 



(56) 



To see how renormalons cancel in physical quantities, let us consider the I/ttlq 
expansion for pseudoscalar and vector meson decay constants in the HQET. It 



reads 31 



fpJm^ 



Ci{mQ/fi) + C2(mQ//i) F(/i) ll + 



Giifi) -b{mQ/fi)- 



+ 



6 



Cmagimq/fi) G2{n) - B{mQ/^) — 



fvVm^ 



rriQ 

Ci(mQ/^)F(/i)|l 
2 



1 



Giifi) -b{mQ/fi) 



A 



+ Oil/ml) 



Crm^gimq/fi) G2(/i) - B{mQ/^) — 



+ Oil/ml), (57) 



where Cmag("^Q/At) = l + 0{as) is the coefficient of the chromo-magnetic operator 
in the effective Lagrangian (p4|), while B{mQ/^) = 1 + 0{as) and b{mQ/fi) = 
0{as) are coefficients that appear at order l/rriQ in the expansion of the currents 
[0. F(/i), G'i(/i), and 6*2 (a'-) are hadronic parameters, which are independent of 
mq. Both G'i(/u) and 6*2(1^) mix with lower dimensional operators and contain 
ultraviolet renormalons, as can be seen from the renormalization-group equations 



d/i 



Cmag(mQ/^) G2(^)] = Y2 ^ d~ ^^'^^/^^ ■ 



(5^ 



Requiring that in (^) the infrared renormalons in the coefficient functions cancel 
against the ultraviolet renormalons in A and Gi{fi), we obtain the relations 



Ci + (72 6mQ mg ^ 



C1 + C2 
ACi 



6?71, 

Ampoio 



Ci GniQ 

To order 1/Nf, they simplify to 



(& - fi) - — (AGi - 2C^ag AG2 



(59) 



AG, = -^ (4ACi + AC2) + 0(iV72) = ^^^ + 0(iV72) , 



AGs = -^!J^ _ !^ AC2 + 0(Ar-2) 



12 



8 



_^:':^Mi_^o{Nf] 



6 



(60) 
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This determines the ultraviolet renormalon ambiguities in the hadronic param- 
eters Gi{fi). The situation encountered here is general for heavy-to-light tran- 
sitions; since there are always at least two hadronic parameters that contain 
ultraviolet renormalons, it is not possible to derive a consistency condition for 
the infrared renormalons in the coefficient functions Ci and C2. However, the 
residues of the renormalons in the coefficient functions determine in a unique 
way the residues of the ultraviolet renormalon poles in the hadronic parameters 
of the HQET. 

5 Inclusive Decays of Heavy Hadrons 

After the analysis of exclusive transitions, we will now consider current-induced 
inclusive decays of hadrons containing a heavy quark. Examples are the semilep- 
tonic decays B ^ Xqii? and Ai, -^ Xqi u, where g = c or m, as well as the rare 
decay B ^ Xs'~f. The flavour- changing current relevant to semileptonic decays is 
g7'^(l — 75) 6. For the penguin-induced transitions, it is of the form gcr'^'^(l±75) b. 
The inclusive decay distributions can be calculated in powers of 1/mi, using an 
OPE for the transition amplitude ||3^-||40|1 



T(v,p) = -i f d^xe-'"'-'' {H{v)\T {b{x)riq(x),q{0)r2b{0)]\H(v)) , (61) 

where H{v) denotes the decaying 6-flavoured hadron with velocity v, p is the 
momentum carried by the current (in the cases above, the total lepton or photon 
momentum, respectively), and Fj are abbreviations for the appropriate Dirac 
matrices. The OPE is constructed by performing a phase redefinition [cf. (|l^)] 

by{x) = exp(— imfcf ■ x) b{x) = h^{x) + 0(l/mf,) (62) 

to pull out the leading dependence of the fields on the heavy quark mass. The 
next step is to write T{v,p) as a sum of coefficient functions multiplying local, 
higher dimensional operators. The coefficients are determined by evaluating the 
diagrams shown in Fig. §, where the momentum of the 6-quark is as usual written 
in the form pb = rrih v + k. The residual momentum k is equivalent to a derivative 
acting on the rescaled heavy quark field b^. 

We will evaluate the contributions in the OPE including terms of order l/rub 
and '^-/Nf. In general, the equation of motion can be used to relate all terms of 
order l/rrib to the residual mass term in the HQET Lagrangian, which is itself of 
order l/Nj. Hence, it will be sufficient to evaluate the l/mi, corrections at tree 
level. Let us then start with the discussion of the tree diagram in Fig. |^. It gives 

^*- = ^ ^1 V w^ — r-^ ^2 ) , (63) 

mbf — p — rUq + lip 
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where we use the short-hand notation {H{v)\byrby\H{v)) = (F). The tree 
diagram contains the propagator of the g-quark in the background field of the 
hght degrees of freedom in the decaying hadron. To proceed, we expand the 
propagator as 

1 1 1^1 

+ ... , 



^^fe f ~ ll ^ ^^q + '^P '^b f '^ 1^ ^ ^ "^fe i ~ li ~ "^12 "^b f ^ 'fl ^ ""^g 

(64) 
where the second term is of order l/rrife relative to the first term. The forward 
matrix element of any local operator h^ T^ iD"" h^ containing a single covariant 
derivative can be evaluated, up to l/ruh corrections, using (|^) together with the 
equation of motion iv ■ D hy = Smhy, where 6m is the residual mass term. It 
follows that 

(r,zD") = 5m(r„i;") + ... , (65) 

where Ta denotes an arbitrary Dirac matrix, and the ellipses represent terms that 
are suppressed by one power of l/nib. Applying this relation, and resumming the 



expanded propagator (Q), we find 



Ttree = ( T, ^ T,) + . . . , (66) 

[rrih + dm) f ~ ff ~ m,q 

where the ellipses represent terms of order l/m,l relative to the leading term. We 
observe that, as in the case of exclusive decays, the residual mass term always 
appears together with the HQET expansion parameter rrih in the combination 
m^° '^ = m^h + Sm., i.e. it is the pole mass that enters the tree-level expression for 
the transition amplitude. The infrared renormalon in the pole mass leads to an 
ambiguity given by 

ATtoee = ( Ti -^- (-Ampoie f^) -^ T^ ) . (67) 

m^bf — p — m,q m^hf — p — m,q 

To see how this renormalon is cancelled, let us now turn to the calculation of 
the radiative corrections depicted in Fig. |^. We study the Borel transform of the 
transition amplitude to order 1/A^/ using the resummed gluon propagator (|13D . 
In the calculation, we only keep terms that have a renormalon pole at u = 1/2. 
We obtain: 

%erte.{u) = ^ e-""/' fiTil - 2u) {Ti —^ ^ —^ T^) + . . . , 

Po mbp-p-mq mbp-p-mq 

Tbox(n) = -^e-^/Vr(i-2w)(ri ^ ^ ^ T2) 



/?o rribf — ^ — m,q WLbf — ^ ~ mi 

+ ^e-^/^Ar(i-2«)(r,^ r,) + ... 

(3q mb mbp-p- ~" 

6Cp 
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rriq 



twFR(n) = --f e-^/2 ^ r(l - 2u) ( V, -— V,) + . . . . (68) 

Po 'Tife mbp -p~mq 
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The ellipses represent terms that are regular at u = 1/2, and terms of order 
1/Nj. Note that there is no renormalon contribution from the renormalization 
of the g-quark propagator. Moreover, the renormalon poles with residues pro- 
portional to the tree diagram cancel between the box graph and wave-function 
renormalization. For the sum of all loop contributions, we find 

TioopsM = ^e-^/Vr(l-2«)(ri— -^- ^—r^ T,) + . . . . 

Po mbf-p-nig mbp-p-niq 

(69) 

From the residue of the pole at u = 1/2, we obtain for the renormalon ambiguity 

ATioops = ( Ti ^ Ampoic ^ -^ Fa ) + 0{Nf) . (70) 

f^bf ~ ff — f^q rribf) — ff — rriq ■' 

As expected, the sum of all contributions in the OPE for the transition amplitude 
is free of renormalon ambiguities: 

AT = ATt,ee + ATioops = . (71) 

That this cancellation occurs was conjectured by Bigi et al. in Ref. 0, however 
without presenting an explicit calculation. In fact, it was claimed that infrared 
renormalons only appear in the vertex corrections and mass renormalization, 
but not in the box diagram. Our calculation shows that this is not correct .0 
Nevertheless, we confirm that the cancellation occurs when all diagrams are taken 
into account. 

The situation encountered here is special in that to order l/mb there do not 
appear non-perturbative corrections when the pole mass is used in the OPE of the 
transition amplitude. Hence, at this order there are no ultraviolet renormalons. 
What we have demonstrated above is a cancellation of infrared renormalons. 
Consider, as an example, the total decay rate for the process B — > X^ i v- It can 
be calculated from the imaginary part of the transition amplitude. Neglecting 
the mass of the w-quark, one obtains the well-known result ETl 



F(5 ^X^lv) = ^['^"f C{m,) {l + 0(l/m^)} , (72) 
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where 
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C(m,) = {mrr (1 - ^^1^ (^TT^ - -j + . . . I . (73) 

We have shown that the infrared renormalon at m = 1/2 in the pole mass is 
cancelled by an infrared renormalon in the perturbative series. It is possible to 
eliminate these renormalons explicitly by introducing a heavy quark mass m^ 
renormalized at short distances instead of using the pole mass 0, ^ . As long as 

^We note that the renormalon contributions in the individual diagrams are the same in all 
covariant gauges. 
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mf differs from m^°^ by a multiplicative factor Z[as(^fe)], this substitution does 
not induce l/rrih corrections to the decay rate ([72|) . For instance, we may work 
in the MS scheme and use the running mass mf,{fi) evaluated at /i = rrif,. This 
leads to 

CW = K(,n.)f{l-?^(.^-f )+.,,}. (74) 

The perturbative series in this expression does no longer contain a renormalon 
at M = 1/2. Note, however, that at some higher order in the 1/mf, expansion 
there will appear ultraviolet renormalons in the non-perturbative corrections to 
the decay rate (|7^ . Correspondingly, the coefficient C{mb) must contain infrared 
renormalons at larger values of the Borel parameter u, which cannot be eliminated 
by introducing the renormalized mass mf-. 

6 Summary and Conclusions 

We have investigated the appearance of renormalons in the HQET by considering 
the l/mq expansion for exclusive heavy-to-heavy and heavy-to-light transitions, 
as well as for inclusive decays of heavy hadrons. We have argued that, in general, 
infrared renormalons in the coefficient functions of HQET operators are compen- 
sated by ultraviolet renormalons in the matrix elements of higher dimensional 
operators, and we have identified which of the HQET matrix elements contain 
such ultraviolet renormalons. In the case of heavy-to-heavy transitions, the sym- 
metries and the equation of motion of the effective theory lead to five consistency 
relations among the infrared renormalons in the pole mass and the coefficient 
functions. We have checked that these relations are satisfied to next-to-leading 
order in an expansion in powers of l/iV/- 

The most important, though not surprising, result of our analysis is that the 
appearance of renormalons does not alter the structure of the heavy quark ex- 
pansion, and does not invalidate any of the predictions derived using the HQET. 
In particular, Luke's theorem, as well as relations between weak decay form fac- 
tors, remain valid. In this sense, there is no "renormalon problem" in the HQET. 
However, as in any OPE it is true that some of the dimensionful hadronic param- 
eters describing the non-perturbative corrections in the heavy quark expansion 
have an intrinsic uncertainty of order Aqq^^. An example is provided by the mass 
parameter A. In the practical form of the OPE, in which dimensional regular- 
ization is employed in the calculation of the coefficient functions, ambiguities 
arise from the necessity to specify a resummation prescription to regulate the 
divergent asymptotic behaviour of perturbation theory. In the literal form of 
Wilson's OPE, they arise from the introduction of a hard factorization scale /i. 
The hadronic parameters of the effective theory then exhibit a power-like depen- 
dence on /i, in a way that depends on how the cutoff is implemented. In both 
cases, to define these parameters precisely would require one to fix terms in the 
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coefficient functions that are exponentially small in the coupling constant. As 
long as one works with truncated perturbative expressions for the Wilson co- 
efficients, the errors due to the truncation are parametrically larger than power 
corrections. This type of ambiguity is inherent in any OPE and as such cannot be 
avoided. In this context, we note that the introduction of a short- distance mass 



instead of the pole mass, which was proposed in Refs. 0, ^ , does in general not 
help to eliminate renormalons. An exception is the case of inclusive decays of 
heavy hadrons, where this procedure eliminates the leading infrared renormalons. 
On the other hand, such a choice of the heavy quark mass destroys the flavour 
symmetry of the effective Lagrangian of the HQET and is thus unattractive, at 
least in processes that involve more than one heavy quark flavour. 

Finally, we like to point out that our somewhat formal investigation of renor- 
malons can serve for tests of HQET calculations. In some cases, the requirement 
that a compensation of infrared and ultraviolet renormalons occurs leads to non- 
trivial relations. An example is provided by the consistency conditions (|36|) for 
heavy-to- heavy transitions. Using a similar argument, we could show with an ex- 



plicit calculation that a sum rule derived by Shifman et al. ||2^, which has been 
used to put a bound on the hadronic form factor that enters the extraction of | Vcb\ 
from semileptonic decays, must be incorrect. This sum rule relates a physical ob- 
servable to a theoretical expression in which infrared renormalons in a coefficient 
function do not match with ultraviolet renormalons in non-perturbative parame- 
ters. This expression has an intrinsic ambiguity and thus cannot be complete. A 
further investigation of what goes wrong with the argument presented in Ref. [^ 
is necessary before any useful phenomenological bound can be derived. 

While this paper was in writing, we became aware of a preprint by Beneke 
et al. 1^, who demonstrate the cancellation of infrared renormalons in inclusive 



decay rates. Their results agree with our Sect. ^ 
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Figures 




Figure 1: Borel transform of the heavy quark self-energy to order 1/Nf. 
The resummed gluon propagator (|1^) is denoted by the dashed bubble. 
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Figure 2: Vertex contribution to the matching calculation of the coeffi- 
cient functions of heavy-heavy currents. 
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Figure 3: Tree- level contribution and radiative corrections to the tran- 
sition amplitude. 
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